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The set of solutions inferred by the generic maximum entropy (MaxEnt) or maximum relative 
entropy (MaxREnt) principles of Jaynes - considered as a function of the moment constraints or their 
conjugate Lagrangian multipliers - is endowed with a Riemannian geometric description, based on 
the second differential tensor of the entropy or its Legendre transform (negative Massieu function). 
The analysis provides a generalised least action bound applicable to all Jaynesian systems, which 
provides a lower bound to the cost (in generic entropy units) of a transition between inferred positions 
along a specified path, at specified rates of change of the control parameters. The analysis therefore 
extends the concepts of "finite time thermodynamics" to the generic Jaynes domain, providing a link 
between purely static (stationary) inferred positions of a system, and dynamic transitions between 
these positions (as a function of time or some other coordinate). If the path is unspecified, the 
analysis gives an absolute lower bound for the cost of the transition, corresponding to the geodesic 
of the Riemannian hypersurface. The analysis is applied to (i) an equilibrium thermodynamic 
system subject to mean internal energy and volume constraints, and (ii) a flow system at steady 
state, subject to constraints on the mean heat, mass and momentum fluxes and chemical reaction 
rates. The first example recovers the minimum entropy cost of a transition between equilibrium 
positions, a widely used result of finite-time thermodynamics. The second example leads to a new 
minimum entropy production principle, for the cost of a transition between steady state positions of 
a flow system. 



1. INTRODUCTION 

Jaynes' maximum entropy principle (MaxEnt) and its extension, the maximum relative entropy principle 
(MaxREnt), based on the principles of inductive (probabilistic) rather than deductive reasoning, arguably consti- 
tutes one of the most important tools for the solution of indeterminate problems of all kinds [TJ [2j |3j HJ [5l El H] • ^ n ^ s 
method, one maximises the entropy function of a system - a measure of its statistical spread over its parameter space 
- subject to the set of constraints on the system, to determine its "least informative" or "most probable" probability 
distribution [I] [21 E] • By a series of generic "Jaynes relations" , this can then be used to calculate the macroscopic 
properties of the system, providing the best (inferred) description of the system, subject to all that is known about 
the system. Since its inception half a century ago, the MaxEnt and MaxREnt principles have been successfully ap- 
plied to the analysis of a diverse range of systems, including in thermodynamics (its first and foremost application), 
solid and fluid mechanics, mathematical biology, transport systems, networks, economic, social and human systems 

[DEJElSlElEJtzlIHlE]. 

The aim of this study is to examine a valuable extension to Jaynes' generic approach, by endowing the set of solutions 
inferred by Jaynes' method - considered as a function of the set of moment constraints and/or their conjugate 
Lagrangian multipliers - with a Riemannian geometric interpretation, using a metric tensor furnished directly by 
Jaynes' method. The analysis leads to a generalised least action bound applicable to all Jaynesian systems, which 
provides a lower bound for the cost (in generic entropy units) of a transition between different inferred positions of 
the system. The analysis therefore extends the concepts of "finite time thermodynamics", developed over the past 
three decades [I(3[lTJ[l2[I2[ll[T3 

[37l [38l [39l [40l [4TI [42l [43] , to the generic Jaynes domain. The analysis reveals a deep, underlying connection between 
the essentially static manifold of stationary positions predicted by Jaynes' method, and lower bounds for the cost of 
dynamic transitions between these positions. 

The manuscript proceeds as follows. In the theoretical principles of Jaynes' MaxEnt and MaxREnt methods 
are discussed, followed by an appraisal of a generalised free energy (generalised potential) concept associated with 
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Jaynes' method. In the concepts of a Riemannian metric, arc length and action sums and integrals are developed 
in a generic Jaynesian context, leading to a generic least action bound for transitions on the manifold of Jaynes 
solutions. Considerations of minimum path lengths, involving calculation of the geodesic in Riemannian space, are 
also discussed. In Q the foregoing principles are applied to (i) an equilibrium thermodynamic system subject to 
mean internal energy and volume constraints, and (ii) a flow system at steady state, su bjec t to constraints on the 
mean heat, mass and momentum fluxes and chemical reaction rates. The first example (j |4.1| ) recovers the minimum 
entropy cost of a transition between equilibrium positions, a widely used result of finite-time thermodynamics. The 



second example (£4.2) leads to a new minimum entropy production principle, for the cost of a transition between 
steady state positions of a flow system. The analyses reveal the tremendous utility of Jaynes' MaxEnt and MinXEnt 
methods augmented by the least action bound, for the analysis of probabilistic systems of all kinds. 

2. JAYNES' GENERIC FORMULATION (MAXRENT) 

2.1. Theoretical Principles 

The usefulness of Jaynes' method for statistical inference arises from its generic formulation, first expounded by 
Jaynes and other workers in the context of information theory [l] [2j [3] HJ (H [7] , but which can be reinterpreted using 
a combinatorial framework (the "Boltzmann principle") [44[ |45] |46j |47[ [48 1 49, 50j[5l]. In consequence, the method 
can be applied to any probabilistic system involving the allocation of entities to categories; this includes - but is not 
restricted to - thermodynamic systems. For maximum generality, it is useful to include source or prior probabilities qi 
associated with each category i = 1, s, to give the maximum relative entropy (MaxREnt) or minimum cross-entropy 
(MinXEnt) principle. In the event of equal ^, this reduces to the special case of Jaynes' maximum entropy (MaxEnt) 
principle [TJ |2J [31 II [51 [7]. 

The MaxREnt method proceeds as follows. To infer the "least informative" or "most probable" distribution of a 
probabilistic system, we wish to identify its observable realization or macrostate of maximum probability P. This 
is equivalent to maximising the following dimensionless function, chosen for several "nice" mathematical properties 

@a[45]: 

f, = llnP, (1) 

For a system of N distinguishable entities allocated to s distinguishable categories, it can be shown that the distribution 
is governed by the multinomial distribution P = TV! YV i=1 q™* /nj, where rti is the occupancy of the zth category and 
N = J2l =1 rii. In the asymptotic limit N — ► oo, reduces to the relative entropy function [2 (the negative of the 
Kullback-Leibler function [52 , 53 ): 

*> = -5>ln^ (2) 

i=i Qi 

where pi = rii/N is the frequency or probability of occupancy of the ith category. Maximisation of ([2| is subject to 
the normalisation constraint and any moment constraints on the system: 

s 

J> = !' ( 3 ) 

i=l 

s 

^2Pifri = (fr): r=l, (4) 
i=l 

where f r i is the value of the property f r in the ith category and (f r ) is the mathematical expectation of f r {. Applying 
Lagrange's method of undetermined multipliers to Q-Q gives the stationary or "most probable" distribution of the 
system (denoted *): 

R R 
-A - E Kfri 1 - E Kfri 

Pi = Qie = —qie 

(5) 

\ ^ - E Kfri 

Z = e Xo = l^qie 
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where A r is the Lagrangian multiplier associated with the rth constraint, Z is the partition function and Ao = In Z is the 
Massieu function [4]. In thermodynamics, the constraints (f r ) are usually taken to represent conserved quantities, and 
thus correspond to extensive variables (e.g. internal energy, volume and numbers of particles), whilst the multipliers 
A r emerge as functions of the intensive variables of the system (e.g. temperature, pressure and chemical potentials). It 
is useful to preserve this distinction between extensive and intensive variables, even beyond a thermodynamic context. 

By subsequent analyses [TJ [2j [3j HJ |5j El [7J [54] , one can derive the maximum relative entropy fj* and the derivatives 
of i}* and Ao for the system: 

R 

S* = A + X)M/r> (6) 

r=l 

dSj* 



d(fr) 

d 2 Sj* d\ r 



K (7) 

Imr G g (8) 



9(fm)d(f r ) d(f m ) ■ 

#A ° - (frM - («</».) = -W- = -W e -1 ( 10 ) 



n\ r-vA \J I J III I \J I I \ J Tib I 

oX m dX r oX ri 



The second derivatives of Ao in ( 10) express the dependence of each constraint on each multiplier, and therefore give the 
"capacities" or "susceptibilities" of the system (e.g. in thermodynamics, they define the heat capacity, compressibility, 
coefficient of thermal expansion and other material properties [I2j [55j |56] ) . Their matrix 7, the variance-covariance 
matrix of the constraints (with change of sign), is equal to the inverse of the matrix g of second derivatives of 5}* in 
([8|, yielding the generic Legendre transformation between the (/b), •••) and Ao(Ai, A2, ...) descriptions of the 

system [2]: 

g7 = I, (11) 



where I is the identity matrix [2]. From (|8j) or (10) and the equality of mixed derivatives, we also obtain the generic 



reciprocal relations d(f r ) /d\ m = d(f m ) /dX r for the system. 

Jaynes also showed that the incremental change in the relative entropy can be expressed as pQ: 

R R 

dtf = X r ( d (fr) ~ (dfr)) = ^ X ^Q r (12) 

r=l r=l 

where 5W r = (df r ) = Yli=iPidfri an d SQ r = J2i=i d Pifn can be identified, respectively, as the increments of 
"generalised work" and "generalised heat" associated with a change in the rth constraint, and <$(•) indicates a path- 
dependent differential. Eq. ( [l2| ) gives a "generalised Clausius equality" [57], applicable to all multinomial systems in 
the asymptotic limit. 

It is again emphasised that the above relations (p5])-([l2|) apply to any probabilistic system of multinomial form, 
in the asymptotic limit. Although originally derived in thermodynamics, the above-mentioned quantities need not 
be interpreted as thermodynamic constructs, but have far broader application. Furthermore, the relations (|5|-(12) 
apply to the stationary position of any multinomial probabilistic system. The derivatives (fr])- (flOj) therefore relate 
to transitions of the system between different stationary positions, or in other words, to paths on the manifold of 
stationary positions. Whilst the lack of inclusion of non-stationary positions may seem unnecessarily restrictive, such 
geometry provides a sufficient foundation for most of engineering and chemical equilibrium thermodynamics. As will 
be shown, it is also useful for the analysis of many other systems of similar probabilistic structure. 

2.2. The Generalised Free Energy Concept 

It is instructive to insert (12) into the differential of ([6| and rearrange in the form: 



dcj) = -dX = -d\llZ = J2 X rSW r + J2 dX r(fr) = ~dS)* + ^ X rd(f r ) + ^ dX r(fr) (13) 
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The negative Massieu function — Ao is therefore equivalent to a potential function <fi which captures all possible changes 
in the system, whether they be in the entropy, constraints or multipliers. For constant multipliers, it simplifies to the 
weighted sum of generalised work on the system. It thus provides a dimensionless analogue of the free energy concept 
used in thermodynamics. For constant multipliers, (j)\{\ r ] also provides a measure of the dimensionless "availability", 
or the available "weighted generalised work", which can be extracted from a system. By extension of the principles 
of equilibrium thermodynamics, we can thus adopt the potential <p as a measure of distance from the stationary 
state. The system will converge towards a position of minimum </>, representing the balance between maximisation of 
entropy within the system 5}*, and maximisation of the entropy generated and exported to the rest of the universe by 
the transfer of generalised heats SQ r (see [63] for further discussion). The advantage of Jaynes' generic formulation 
is that (j) can be defined for any multinomial probabilistic system, and is not restricted to thermodynamic systems 

pas o nii]. 

Returning to the second derivatives in the last section, we see that Ao can be replaced by —(f) in (|9|)-([To|). The 
latter provides a clean (non-negative) Legendre transformation between matrices g and 7, and thus between the 
(/b)? •••) and 0(Ai, A2, ...) representations of a system. 



3. RIEMANNIAN GEOMETRIC CONCEPTS 



3.1. Generalised Riemannian Metrics and Arc Lengths 



Since the time of Gibbs [58] [59] , examination of the geometry of the manifold of stationary positions has been 
of tremendous interest to scientists and engineers. In thermodynamics, this has typically involved analysis of the 
concave hypersurface defined by the Euler relation S f (Xi,X 2 , ...), where S is the thermodynamic entropy and X r are 
the extensive variables, or alternatively of its Legendre transform, the convex hypersurface ^(Yi, I2, ...) or F(Y±, Y2, ■•■)» 
where ip = F/T is a Planck potential function, F is a free energy and Y r are the intensive variables [55] [56]. Such 
interpretations have led to major advances in the understanding and analysis of thermodynamic processes and cycles 
[58] [59]. However, adoption of the Jaynes MaxEnt framework (^2| permits a rather different insight, based on a 
Riemannian geometric interpretation. As will be evident from the previous discussion, this interpretation extends 
well beyond "mere" thermodynamics, forming a natural adjunct of Jaynes' generic formulation (^2|. 

Consider the i?-dimensional hypersurface parameterised by the constraints {(f r )} or their conjugate Lagrangian 
multipliers {A r }, representing the hypersurface of stationary states within the (R + l)-dimensional space given by 
{(fr)}) or (05 {-V})- If the R parameters are linearly independent, the matrices of the second derivatives g ^ 
or 7 (flOj) are positive definite (i.e. x T gx > or x T ^yx > for any non-zero vector x [60] [6TJ [62]). The matrices 
g or 7 can therefore be adopted as Riemannian metric tensors associated with the stationary state hypersurface 
defined by {(f r )} or {A r }, and used to interpret its geometric properties. Indeed, even if the R parameters are 
not always independent, whence g or 7 are positive semidefinite (i.e. x T gx > or x T jx > for x ^ 0), the 
latter can still be adopted as pseudo-Riemannian metric tensors on the stationary hypersurface. This representation 
was first proposed by Weinhold [TO] [TTJ [12] [13] [14], and its implications in terms of a least action bound were 
subsequently developed, largely within a thermodynamic context, by Salamon, Berry, Andresen, Nulton and co- 
workers US HSl 113 HH1 QS1 1213 I2S 1221 1231 I2H I2H I2Z1 I2H1 1221 1301 I2H 1321 1331 I3H and also by Beretta 
[37] l38| [39] , Diosi and co-workers [40] , Crooks and Feng [41] [42] and Brody and Hook [43] . Some theoretical aspects of 
the adopted Riemannian formulation are discussed in Appendix [A] It must be noted that the Riemannian formulation 
replaces - it cannot be used in conjunction with - the traditional convex or concave hypersurface interpretation 
normally used in thermodynamics and information theory [30] . 

Firstly, the Riemannian geometric interpretation provides an intrinsic differential or line element (its square, a 
metric) with which to measure distances along a specified path on the manifold [60, 62] 1 : 




, d (fm)gmrd(f r ) = ^df T gdf, (14) 

\ m,r=l 



, ^2 dXmjmrdXr = y dA T jdA. (15) 

\ m,r—l 



Strictly, this line element is not a first fundamental form in Riemannian geometry 60 62 ; its use as a distance measure is discussed in 
Appendix [A] 
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where / = (/b), (fs)] T and A = [Ai, A2, A S ] T . Integration between points a and b along a path on the 

manifold, defined by the set of increments df or dA, gives the arc length along that path between those points 

H3HT1ISBIB2]: 



J ds% 

a 

b 

Jds, = J 



R 



, \ E d (f> 

a M m ' r=1 



m/ i/rnr 



df T gdf, 



a dA m 7 mr dA r 

\ m,r=l 





/ 



dA 7dA. 



(16) 



(17) 



The shortest such path is known as the geodesic. An infinite number of other paths on the manifold are also possible, 
of longer arc length, as also given by (16) or dl7|). If the manifold is parameterised by some parameter £ - which can, 



(18) 



but need not, correspond to time t - the arc lengths can be expressed in continuous form as: 



L 



E 



J \ 1 d z 

N m > r=1 



d(fm) djJJ 



d£ : 



• T 



E 



dA^ 



dA r 



A 7 Ad£ 



(19) 



where the overdot indicates differentiation with respect to £. 

The symmetry of the Legendre transformation (11) also permits a further insight. From ([8| and ( [To| ), the metrics 
g mr or 7 mr within the intrinsic differentials (14)-(15) can be substituted respectively by dX r /d(f rn ) or d(f r )/dX m , to 
give: 



ds$ 



ds d 



, E d </» 

\ m,r=l 



dX r 



d{fr 



\ 



^dA r d(/ r ) = v / dA^d/ 



r=l 



a ^ ] dA m — dA r 

\ m,r=l 



. ^2dX r d(f r ) = ^dA-df 

\ r=l 



In consequence, the intrinsic differentials are equal, ds = ds^ 

e 



ds^, and so too are the arc lengths: 



/ 



A • / d£, 



(20) 



(21) 



(22) 



From a Riemannian geometric perspective, it therefore does not matter whether one examines a system using its 
(/2)j •••) or A2, ...) representation. The above identities - touched on by several workers 1271 134 ] [40 ) [4Tj 
- are not surprising, since the Legendre transforms S)* and <p both have the character of entropy-related quantities, 
respectively indicating the (generic) entropy of a system and the capacity of a system to generate (generic) entropy 
[63] . The quantity dA • df therefore expresses the second differential of generic entropy produced due to incremental 
changes in A and / (a generalised force-displacement or fluctuation-response relation). For all changes, dA • df > 
must be valid, to preserve a positive definite metric (whence A • / > 0) [12]; this is in sympathy with a generalised 
form of the second law of thermodynamics, namely "each net mean increment of (generic) entropy produced along a 
path must be positive" . 

One further consideration arises from the recognition that most probabilistic systems involve quantised phenomena, 
which can only be approximated by the above continuous representation. For a system capable only of discrete 
increments in line elements As^* or As^ associated with a minimum dissipation parameter A£ (e.g. a minimum 
dissipation time if £ = i), the arc lengths are more appropriately given as [25] : 



A I 



A I 



A I 



= e a^*,, = E v Af v T go A fv = E yfv gf h 



v=l 
M 



v=l 
M 



v=l 
M 



L (f) = E As <f>,v = E V AA ^ T AA„ = E V ^ Iv&v A£ 



(23) 
(24) 




Environment 



FIG. 1: Illustration of Riemannian geometry concepts, for a two- constraint system represented by i3*((/m), (fr)) (for conve- 
nience the environment is shown as horizontal). 



where v is the index of each increment. The last terms in (23)-(24) invoke the finite difference forms f v = Af v /A£ v 
,/A£y, strictly valid only in the limits A£ v — > 0. The two discrete length scales (23)-(24) are again 



or A r 



AA,. 



equivalent, but there will most likely be some discrepancy between their values due to their finite difference formulation. 



3.2. Generalised Action Concepts and Least Action Bound 

A Riemannian geometry can also be examined from a different perspective [TSJ [T71 [25j |33l [41] , discussed with 
reference to Figure[l] the following analysis largely follows [25 j, converted into generic form. Although applied to 5}*, an 
analogous derivation can be given for the <p representation. Consider a system on the manifold of stationary positions, 
subject to displacements {A(/ r )} in its stationary position. The modified (generic) entropy $)*({(f r ) + A(/r)}) of the 
system can be expanded in a Taylor series about <£)*({ (/r)}) : 



«*({</ P )+A</ P )}) = + £ A r | {</r)} A(/ r ) + i £ dl f m )Q {fr) 



r=l 



A(/ m )A(/ r ) 



{{fr)} 



<9 3 £* 



1 R 

3! m 5 =1 d(ft)d(f m )d(f r ) 



(25) 



A(/,)A(/ m )A(/ r ) 



{</r» 



where use is made of Q. The corresponding change in entropy of the "reservoir" or "environment" of constant 
{\f. nv }, by which this change is effected, is given (exactly) by [17} 123]: 



^({(/ r )+A(/ r )})=^({(/ r )}) + £An { < /r>} A(/ r ) 



(26) 



At the stationary state, A r = X^ nv , whilst from the constraints (conservation laws), A(/ r ) = —A{f r ) env [25_. Addition 
of (25)-(26) thus yields the total change in the entropy of the system and environment for the step process: 



AS)* 



1 R 



d 2 Sj* 



2! 



% d(f m )d(f r ) 



R 

A(/ m )A(/ r ) + - 



3k* 



d s f> 



{</.)} 



f^ =1 d(fe)d(f m )d(f r ) 



A(f e )A(f m )A(f r ) 



(27) 



{</>•>} 



7 



Provided the manifold is smooth, continuous, continuously differentiate (i.e., there are no phase changes in the 
neighbourhood) and the step sizes {A(/ r )} are small, we can neglect the higher order terms in (27), giving: 



Afiv ~ ^ E A </™>- 9mr,v\{Ur)} A (fr)v = J Sv Af v (28) 



m,r=l 



where the subscript denotes the vth equilibration step. In the <\> representation, the analogous form is obtained (in 
this case, giving the loss in </>): 



1 R 1 

A&v ~ 2 ^2 AA ™>^ lmr,v\{K] A Kv = ~ AA V T J v AA V (29) 



m,r=l 



The total increase in entropy or decrease in potential of the system and environment subject to a M-step process is 
therefore: 



M M . M 



Atf tot = f2^K~J2 \ A fv T Sv A/. = \fJ Sv f v (30) 

V = 1 V = 1 'L> = 1 

M M I M 1 T 

- A(/w = - J] A ^ ~ J2 o AA ^ T ^ AA ^ = 9 ^ ^ A & A ^ ( 31 ) 



2 " ^— / 2 



Recognising Af^ as the minimum dissipation parameter for the ^th step (e.g. the minimum dissipation time if £ = t), 
one such term may be factored out of each sum, to give the mean minimum dissipation parameter e n for n G {i}*, <f\ 
HUES]. This gives A£* ot = e^*^* or -A0 tot = e J" , with: 



M 1, T 



^*=E^ g-/. A ^ (32) 



M 



v=l 

The summands i f v g v f v or | ^ v A v in d32b - d33b can be viewed as generalised energy terms, akin to the kinetic 
energy in mechanics, with the metric or 7,, representing the "mass" and f v or the "velocity" [64]. The terms 
J ^ or can then be interpreted as the discrete generalised action of the specified process [41 , again by analogy 
with mechanics 2 . From the previous considerations (^3|, the two action sums are equivalent, although once again, 
discrepancies may emerge from their finite difference formulation. 
From the discrete form of the Cauchy-Schwarz inequality: 

(eL- j )(e1^)>-(e1^) 2 (-> 



with a v = \l f T g v f A^ v or y A T j v A A£ v and b v = 1, it can be shown that 



enJn > ^ (35) 



Physically, the number of steps is equal to M = ^max/tn, whence (35) reduces to [T7] [25]: 

or J n > — 

* <>max * S' 



e„Jn>^^- or J n >-^- (36) 



2 Crooks |41| applies the terms "energy" and "action" interchangeably; we consider that the present definitions are more in keeping with 
those used in mechanics. Many authors include the e n term within J7" n , but we here wish to preserve the mathematical structure of a 
generalised action principle. 
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Eqs. (35)- (36) can be considered a generalised least action bound [41 j, applicable to all probabilistic systems amenable 
to analysis by Jaynes' method. Its physical interpretation is that it specifies the minimum cost or penalty, in units of 
dimensionless entropy per unit £, to move the system from one stationary position (at £ = 0) to another (at £ = £ ma ^) 
along the given path at the specified rates A and/or /. If the latter rates proceed infinitely slowly, the lower bound 
of the action is zero, indicating that the process can be conducted at zero cost; otherwise, it is necessary to "do 
generalised work" to move the system along the manifold of stationary positions within a finite parameter duration 

£>max • 

The generalised least action bound thus provides a lower bound for the "transition cost" of a process (in entropy- 
related units). If the process is reversible, the cost would be zero, but no process can be reversible in practice. 
Identification of this minimum cost is of paramount importance: there is no point in undertaking expensive changes to 
the process, or initiating costly social or political changes, in the attempt to do better than the minimum predicted by 
(35)- (36). Taking a thermodynamic example, the method can be applied to determine the minimum cost of industrial 



processes such as work extraction from combustion, a question of fundamental importance to human society. Most 
thermodynamics and engineering textbooks give the Carnot limit as the theoretical limit of efficiency, but the limits 
imposed by finite time thermodynamics are more restrictive (see £4.1). 

The generalised least action bound therefore emerges from the Riemannian geometry of the state space, and hence 
from somewhat different considerations than the principle of least action employed in mechanics [63] . We consider 
that the two principles are connected, but are unable to examine this topic further here. For further exploratory 
expositions, the reader is referred to the work of Crooks [41], Caticha [65] and Wang [67] [68]. 



The above discrete sums (30)-(31) can also be presented in integral form. Consider a system represented by f)*, 
subjected to a finite change in the multipliers AA r due to movement of the reference environment. The incremental 



change in entropy is, again to first order (compare (27)) [17, 25, 33 , [4Tj : 



1 R 



(37) 



R 



Substituting AA r = ^ 9mr^(fm) from (M), and assuming a first order decay process: 



{fm) 



{fm) {fm)env ^{fm) 



where e^* is a minimum dissipation parameter (reciprocal rate constant), (37) yields: 



1 R 

dS)* = - ^ (fm)9mrd(f r )ex 



(38) 



(39) 



The total change in entropy Afi* ot = j^ max d$)* is then obtained as: 

^max ^max 



(40) 



Similarly, in the (j) representation, we obtain: 



\ A T 7Ae di 



(41) 



In the continuous representation, the process does not proceed by a series of finite steps; instead, the reference 
variables continuously move ahead of those of the system [17] [25j [33]. However, we still see the influence of a finite 
decay parameter e n , which on integration yields the mean parameter e n . Each J n term above can be regarded as the 
action integral corresponding respectively to (32)- (33). Based on the integral form of the Cauchy-Schwarz inequality 
[62] , it can be shown that the integral actions also satisfy the least action bound (35)- (36), with L n in integral form 
[171 [25] [33] [41]. 



Finally, for the least action bound (36) to achieve equality, the summands or integrands of the arc length L n and 



action J n must be constant. This gives the simple result that for slow processes with constant dissipation parameter 
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e n = e n , the minimum action (whence minimum in e n J n ) is attained by a process which proceeds at a constant speed 

$n = \[~f g / = \l A [28j [23 [33]. For a constant metric, this is equivalent to constant rates of change of the 
parameter vector / and/or A. For systems with a variable dissipation parameter e(£), it was first considered that the 
minimum is attained at the constant speed ds/dn, expressed in the "natural" parameter units n = £/e [28 , 29 , 33, 40J. 
This however oversimplifies the minimisation problem, which is better handled within a discrete (stepwise) framework 
[25j[35]. As discussed in £4.1, such principles have been widely applied to thermodynamic systems. 



3.3. Minimum Path Length Principle 



The above discrete or continuous forms of the least action bound (35)-(36) are based on consideration of a specified 



path on the manifold of stationary positions, of arc length L n . In many situations, we may wish to determine the 
path of minimum arc length L n ^ min - the geodesic - on the manifold of stationary positions. From the calculus of 
variations, this is given by the Euler-Lagrange equations [66] : 



ds 



d ds 



df 

ds^ 
dA 



d£ df 

d ds^ 







(42) 
(43) 



where s^* = y / 
parameters /, A E 



g / and s 

7J3>2 



A 7 A are the integrands respectively of L^* or ( 16 )-( 19 ) . For two-dimensional 



(42)- (43) can be reduced further in terms of the three unit normals to the surface, giving the 
curve(s) on the manifold for which the geodesic curvature vanishes [43, 60, 62 . Depending on the specified problem, 
a geodesic may not exist, or there may be multiple or il-defined solutions. Provided it does exist, a geodesic leads to 
the double minimisation principle: 



Jn > 



2£« 



> 



L r , 



2& 



(44) 



where the right hand side indicates the absolute lower bound for the action, irrespective of path. This principle has 



been applied to thermodynamic systems, as will be discussed in £4.1 



4. APPLICATIONS 



As noted, the foregoing Riemannian geometric interpretation (§3]) has mainly been presented within an equilibrium 
thermodynamics context [13 [T^ [12 [13 113 12B 

S3 [2J [42j [43] , although it has been applied to non-equilibrium thermodynamic and flow systems [69, 70] I7T | [72 ] [73]. 
information coding [74] and in economics [75]. In the following sections, the utility of Riemannian geometric properties 
and the least action bound are demonstrated for two types of system: a thermodynamic system at equilibrium, and 
a flow system at steady state. 



4.1. Equilibrium Thermodynamic Systems 

The application of Riemannian geometric principles to equilibrium thermodynamic systems has constituted a major 
new development over the past three decades, forming an important plank of finite-parameter or (with £ = t) finite- 
time thermodynamics [2TJ [34] [35] . Such analyses have progressed in four overlapping stages: 

• The initial studies by Weinhold p~3 [11] [12] [T3 E] and early work by Salamon, Andresen, Berry, Nulton and cowork- 
ers [13 [12 [IB > 20, 25, [23 [29] all examined a manifold based on an internal energy representation U(Xi,X2, -•■), as 
a function of extensive variables X r , which include the thermodynamic entropy S. The resulting quantity e~uJu ( m 



the present notation) was interpreted as an availability or exergy function, with (36) indicating the most efficient 
path (defined by the minimum amount of work or minimum loss of availability) required to move the equilibrium 
position of the system [12 HE] . Such analyses complement the thermodynamic geometry used by Gibbs [53 [59] , 
and fit well with the traditional heat- work framework of 19th century thermodynamics. 
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Subsequently, following earlier pioneering works [69l [76] , the entropy manifold 5 , (X L ,X 2 , ...) was examined from a 
Riemannian perspective [20] [2TJ [23] [25] [3TJ [32 1 33 , 33 , 34j [4TJ |42] , where X r are the new extensive variables, of course 
related to the U(X\, X^ ...) representation by Jacobian transformation [20]. The quantity e~sJs was interpreted 
as a measure of energy dissipation or entropy production, again providing a measure of process efficiency. It was 



realised that the lower bound in (36) provides a formal, mathematical definition of the degree of irreversibility 
of a transition between equilibrium positions, with reversibility only for Js = (a definition vastly preferable to 
the cumbersome word-play still used in thermodynamics references; see the scathing criticism by Truesdell [77]). 
However, the primacy of the entropy representation over that based on internal energy was not fully appreciated in 
these early studies. The applicability of Riemannian geometry in other contexts - based directly on the MaxEnt 
framework of Jaynes [1, 2 - is hinted at by Levine [27] , but unfortunately was not developed further at the time, 
nor, to the authors' knowledge, in any subsequent studies. 

» Several studies have considered an entropy representation based on a metric defined on a probability space {pi}, 
either from the Boltzmann principle [76] or using a Shannon or relative entropy measure [22] [23] [24[ [26] [27j |37j [38] 
[39] [HJ E2J |43] . Several authors [37, 38, 39 , [41] [42] 03]) have extended this analysis, to establish a connection with 
the Fisher information matrix [78] and an "entropy differential metric" of Rao [79] . The analysis is also intimately 
connected with paths in a space of square root probabilities, and thence to formulations of quantum mechanics 
[37] [38] [39] . These insights - not examined further here - demand further detailed attention; they may well furnish 
an explanation for the utility of extremisation methods based on the Fisher information function in many physical 
problems [80] . 

» Finally, several workers realised that Riemannian geometric principles can be applied to Legendre-transformed 
representations, e.g. based on various forms of the free energy F (conjugate to U) or the negative Planck potential 
F/T (conjugate to S), as functions of the intensive variables (or functions thereof) [23, 27j[4TJ[42]. This approach 
offers particular advantages for the analysis of real thermodynamic systems, in which the control parameters tend 
to be intensive rather than extensive variables (the canonical ensemble), and for which the intensive variables do 
not exhibit sharp transitions or singularities associated with phase changes, as is the case for extensive variables 



[4Tj . Furthermore, the resulting metric is equivalent to the variance-covariance matrix of the constraints (10), and 
is therefore connected to fluctuation-dissipation processes within the system. 

For completeness, we demonstrate - for a microcanonical thermodynamic system - how Riemannian geometric 
properties emerge as an inherent feature of Jaynes' MaxEnt formulation. Consider an isolated thermodynamic system, 
containing molecules of possible energy levels 6i and volume elements Vj, subject to constraints on the mean energy 
(U) and mean volume (V). We consider the joint probability Pij of a particle simultaneously occupying an energy 
level and volume element, giving the entropy function: 

fieq = ~ ^2^2 Pij Inpzj, ( 45 ) 

i 3 

where, without knowledge of any additional influences, we assume that each joint level ij is equally probable (hence 



the priors cancel out). Eq. (45) is maximised subject to the constraints: 

i 3 

EE^ = <£0, ( 4? ) 

i 3 

EE^ = w, ( 48 ) 

i 3 

to give the equilibrium position: 

~ — XuUi — \ v Vj i 

1 -\uU % -X v V 3 



i 3 



-X u U l -X v V J 



e -AuU z -A v V^ ( 4 g) 



where Z is the partition function. From the existing body of thermodynamics, we can identify the Lagrangian 
multipliers as \jj = 1/kT and Ay = P/kT, where k is the Boltzmann constant, T is absolute temperature and P is 
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absolute pressure. Eq. (49) and Jaynes' relations (|6|-(11) and (13) then reduce to: 

n *. ^Lp-Ui/kT-PVj/kT 



S* = kf>* = k\nZ 



(U) P(V) 
T T 



kA 



eg 



' dS* dS* ' 


T 


"1 p~ 


T 


_d(uy o(v)_ 




T t _ 





d 2 S* 



d 2 S* 

«9(rj) 2 ' d(U)d(V) 
d 2 S* d^S* 



d(V)d(U)' d(V) 2 
iP = k<p eq = -k\nZ = -S* 
dtp dtp 



d 



d , 


( P 




\T 


d , 


(P 


d(V)\ 




_ G 




" T 





' eg 



leg 

k 



d 2 ip 



d 2 tp 

d 2 tp d 2 ip 



.d&aay d(%) 2 



(U) , P(V) 
T T 

= [(U),(V)] T 

d(U) d(V) 

WY W) 

d(U) d(V) 

[dify d($) 



kgeq 



k 



where S* is the thermodynamic entropy at an equilibrium position, ip is the negative Planck potential 
Massieu function [83 ) and G is the Gibbs free energy. By Jacobian transformation of variables, using 
material properties (susceptibilities) [12, 55l 156]: 

Heat capacity at constant pressure: Cp = y ^ J 

Isothermal compressibility: kt 

Coefficient of thermal expansion: 



i fd(vy 

(V) \dP y T 

(V)\ dT ) p 



where (H) = (U) + P{V) is the enthalpy, as well as the equality of cross-derivatives (Maxwell relation): 

d(V) _ d(U) 



the ip metric (56) reduces to' 



.3. 



leg 

k 



T(V) 



whence from (57): 



kgeq = 



-k t P 2 + 2aPT - ^pf, 
k t P - aT, 



T 2 (a 2 T(V) - ktC'p) 



ktP — aT 
—kt 



k t P - aT 



k t P - aT, k t P 2 - 2aPT - 



C P T 
(V) 



(50) 
(51) 
(52) 

(53) 

(54) 
(55) 

(56) 

(57) 

(negative 
the following 

(58) 
(59) 
(60) 

(61) 
(62) 
(63) 



3 The first variance is given erroneously by Callen |55l I56| . 
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Using (52), (55) and (62)-(63), the (dimensional) arc lengths (18)-(19) and action integrals (40)-(41) are obtained as: 



L, 



Js> 



feg kg eq f eq d£ 



2 feq k Seq f pq d£ = 




-Cpt 2 + 2a(V)TTP - k t (V)TP 2 



k 2 T 2 



CpT 2 + 2a(V)TTP - k t (V)TP< 
2k 2 T 2 



d£ 



(64) 
(65) 



-I 



I 



^eq j A- e q d^ 



{V)(U} [k t {U) + 2{V){k t P - aT)] + {V} 2 [P{V){k t P - 2aT) + C P T 



(66) 



T 2 {V){-C P K T + a 2 T{V)) 



1 A 1«i 



A eq d£, 



(V)(U) [k t {U) + 2{V)( KT P - aT)] + (V) 2 [P(V)(k t P - 2aT) + C P T] 
2T 2 (V)(-C p k t + a 2 T(V)) 



(67) 



Using (58)-(60), these two sets of measures can be shown to be equivalent. The above equations must be integrated 
along the particular thermodynamic path followed by the process, as defined by the velocities {T, P} or {(?/), (V')}- 
For a process which follows a pre-determined path, e.g. an adiabatic, isothermal, isovolumetric or isopiezometric curve, 
this can be simplified by expressing the velocities (e.g. P) as functions of one independent velocity (e.g. T). 

To comment on units: if the above quantities were calculated using the "pure" metrics g eq or y eq , in either case 
the line element ds ni arc length L n and the term e n J n would be dimensionless (whence the action is in reciprocal £ 
units). Use of the "natural" metric kg eq , as conducted here, gives the line element and arc length in V JK~ X and 
the action in Ji^ _1 £ _1 , consistent with AS^ ot = ts*Js* being in entropy units. In contrast, use of the "natural" 
metric J e q/k gives the line element and arc length in \/KJ~ l and the action in i^J _1 £ _1 . The latter case can be 



rescued by use of a modified line element ds^ = y kA eq ^ 

_1 and the action in JK -1 ^ 



kA eq d£ - as suggested by (52) and (56) - giving the line 



element and arc length in V JK~ X and the action in JK x £ 1 . Thus in both the S* and ip representations, the least 
action bound ( 36 ) can be used to determine the minimum entropy cost of a transition from one equilibrium position 



to another, along a specified path on the manifold of equilibrium positions. As noted earlier, for slow processes and 
constant e n , this is attained by a process which proceeds at a constant thermodynamic speed s [28J [29j [33j 00] (a 
more general result is available for rapid processes [84]). For variable e n and/or for stepwise phenomena, the process 
should be divided into individual steps placed at equal distances along the arc length traversed by the process, giving 
the so-called "equal thermodynamics distance" principle [34[ [35j [36] . Such considerations have been applied to the 
optimisation of a wide variety of engineering and industrial batch and flow processes, including engine cycles, heat 
engines and pumps, chemical reactors, distillation towers and many other systems. 

A final important point is that the minimum path length (double minimisation) principle ( |44| ) - involving calculation 
of the geodesic - has been applied to the analysis of equilibrium systems. In early work, this bound was established 
by applying the calculus of variations directly to particular thermodynamic problems, without use of a metric [16] [19]. 
More recently, such lower bounds have been examined for particular thermodynamic systems [34] [43] 1100] . In either 
case, for the entropy representation, this method yields an absolute minimum entropy cost AS* > AS^ n for a 
transition between two equilibrium positions at particular rates of change, irrespective of the path. For cyclic or flow 
processes, this therefore gives a minimum entropy production principle S > S m i n , providing one of the key concepts 
of finite-time (or finite-parameter) thermodynamics. 



4.2. Flow Systems 



We now consider a flow system consisting of a control volume, subject to continuous flows of heat, particles 
and momentum, and within which chemical reactions may take place. A few workers have examined such non- 
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equilibrium systems previously within a Riemannian context, including for the Onsager linear regime [69] [70] and 
for extended irreversible thermodynamics [7TJ [72] [73]. A different perspective is provided here, based on a recent 
analysis of a flow system from a Jaynesian perspective [63] . This involves a probabilistic analysis of each infinitesimal 
element of the control volume, which experiences instantaneous values of the heat flux jg j, mass fluxes 2m c °f 

each species c, stress tensor r j and molar rate per unit volume £c d of each chemical reaction <i, where the indices 
X, J ', Cd,N c G {0, ±1, ±2, ...}. We therefore consider the joint probability 7rj = ^x,j,{c d },{N c ] °f instantaneous 
fluxes through the element and instantaneous reactions within the element, giving the (dimensionless) "flux entropy" 
function: 



fist = - ^ 7TJ In TTj , 



(68) 



Again assuming that each joint level / is equally probable, (68) is maximised subject to constraints on the mean 

values of the heat flux (jg), mass fluxes (j c ), stress tensor (r) and molar reaction rates through or within the 
element, as well as by the natural constraint (J3|. This gives the steady state position of the system: 



exp -Cq ■ j QtX 



(69) 



where £q,C c >Ct aim C<2 are the Lagrangian multipliers associated with the heat, particle, momentum and chemical 
reaction constraints, and Z = is the partition function. By a traditional control volume analysis [S5JIB6, 87, 88J, 
the multipliers can be identified as [63_: 



k 

ev 

k 

OVAd 
k T 



T 



F 



(70) 
(71) 

(72) 
(73) 



where \i c is the chemical potential of the cth constituent, M c is the molar mass of the cth constituent, F c is the 
specific body force on species c, v is the mass- average velocity, Ad is the chemical affinity of the dth reaction (< for 
a spontaneous reaction), V is the Cartesian gradient operator, and 6 and V respectively are characteristic time and 
volume scales of the system. Generalising each component of the above multipliers as ( r and constraints as (j r ) with 
r G {1, Jaynes' relations (|6|)-([TT|) and (13) reduce to: 



Sj* st =lnZ + J2( r (j r ) 



9 2 m t 

9(h) 2 



dVi)'""d<j R ) 



n T 



d 2 Sj* st 
d{j R )d{h) 



d(ji)d(j R ) 



d(j R ) 2 
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dCi d( E 
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r aci 
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d'd 
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(74) 
(75) 

(76) 

(77) 
(78) 
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= I 







(79) 



(80) 



where & can be identified as the local entropy production per unit volume (units of JK 



m 



3 s x ). A flow system 



subject to constant flux and reaction rate constraints will therefore converge to a steady state position defined by a 
maximum in the flux entropy and a minimum in the flux potential <j) st . If these effects occur simultaneously, the 
system will converge to a position of maximum <r, therefore providing a conditional, local derivation of the maximum 
entropy production (MEP) principle [63], which has been applied as a discriminator to determine the steady state of 
many non-linear flow systems [55| 150 1 151} 1551 155| 157] ■ 

In Onsager's analysis of transport phenomena in the vicinity of equilibrium [98] [99] , the fluxes and reaction rates 
are considered to be linear functions of the "forces" (the driving gradients and chemical affinities). In the present 
terminology, this would be written as: 



0V) = Kj2L° rm Cn 



(81) 



where l? rra are the (constant) phenomenological coefficients at the zero-gradient position (i.e., at equilibrium) and 
K = k/OV. In the present analysis, we do not claim linearity between (j r ) and £ m , nor consider that the system 
is "close to equilibrium", but simply adopt the partial derivatives d(j r )/d( m within the metric *y st (79) as a set of 
parameters (functions of £ m ) with which to analyse the system. The present analysis therefore encompasses, but is 
not restricted to, Onsager's linear regime. The diagonal and many off-diagonal terms can readily be identified as 
functions of the conductivities (transport coefficients) and chemical reaction rate coefflcents [88 : 



Heat conductivity: k %3 
Diffusion coefficient, species c: D c %3 
Viscosity coefficient: 



d(jci) 

(d<X 

d(T tJ ) 



d 



d 



dt 



Rate coefficient, reaction d: 



dC r . 



dC r 



(82) 
(83) 
(84) 

(85) 



where C c is the concentration of species c (units of kg m 3 ; often used as a proxy for the chemical potential /i c ), 

(C c d) i s the mean rate of change of concentration of species c in the dth reaction (units of kg m -3 s _1 ), v c & is the 
stoichiometric coefficient of species c in the dth reaction (positive if a product), and the indices z, j, i G {x, y, z}. The 
remaining off-diagonal terms consist of the cross-process conductivity coupling coefficients and conductivity-reaction 
rate coefficients. The Riemannian metric j st can therefore be regarded as a function of the material properties or 
susceptibilities of a flow and chemical reactive system, in the same way that the Riemannian metric for an equilibrium 
system j eq is a function of its various susceptibilities, such as Cp, kt and a (£4.1). As with equilibrium systems, 
an abrupt change in a given component j s t,rm with ( m can be interpreted as the boundary of a phase change in the 
system. Notice also that symmetry of ^ st yields a set of Maxwell-like relations [63] : 



9(jr) d(jm) 



(86) 



These apply to all infinitesimal volume elements of a flow system, not merely those in the vicinity of equilibrium. 



Eqs. (86) considerably simplify the set of parameters needed for analysis, from R 2 to (^J 1 ) coefficients; further 



simplifications may be attainable in certain systems due to geometric and tensor symmetries [88] . 
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The above relations (74)-(79) can now be applied to develop a Riemannian description of a flow system on the 



manifold of steady state positions. In terms of the generalised derivatives, the (dimensionless) arc lengths (18)-(19) 
and action integrals (40)-(41) are obtained as: 





£,max £,ma 

J y/A a t'f at d£= J 



A st j st A st d£ 



o 

Smax 



(87) 



&Cr d(j r 



Jst — 







where, as shown, the two alternative $)* st and (j) st measures are equivalent. Once again, these equations must be 
integrated along the particular path taken between the initial and final steady state positions. 

To comment on units: since the above quantities are calculated using the "pure" metrics g st or y st , the resulting line 
element ds st , arc length L st and the term e st J s t are dimensionless. Use of the "natural" metric Kg st , for K = k/OV, 
therefore gives the line element and arc length in V JK~ 1 m~ 3 s~ 1 and the action in JK~ 1 m~ 3 s~ 1 £ ) ~ 1 , thereby giving 
tstJst m units of entropy production per unit volume. Similarly, use of the "natural" metric ^/ st / K in conjunction 
with the dimensional constraint vector KA st gives the line element and arc length in V JK~ 1 m~ ?J s~ 1 and action 
in JK~ 1 m~ 3 s~ 1 £ > ~ 1 , again giving e st J s t in units of entropy production per unit volume. The least action bound 



( 36 ) therefore yields a minimum entropy production principle, which sets a lower bound for the entropy production 
associated with movement of a flow system from one steady state position to another along a specified path. From 
the previous analysis, this involves two separate minimisation principles: 

• If the path is specified, the process of minimum entropy production will be one which proceeds at constant speed s, 
assuming a slow process and a constant dissipation parameter e. Alternately, if the dissipation parameter e is not 
constant, the minimum entropy production process will be given by a constant arc length speed, in accordance with 
a steady state analogue of the "equal thermodynamic distance" principle [25j [34] [35] [36] . 

• If the path is not specified or can be varied, an absolute lower bound for the entropy production is given by the 



geodesic in steady state parameter space, in accordance with the methods of £3.3 



Although they share a similar name, the minimum entropy production principle derived herein is quite different to 
that of Prigogine [86], which concerns the selection of a steady state position relative to possible non-steady state 
positions, and which only applies to the Onsager linear regime. Similarly, it differs from the minimum entropy 
production principle obtained by t he ap plication of Riemannian geodesic calculations to the manifold of equilibrium 
positions, discussed at the end of £4.1 [T6l [19j [34j [43j 1100] . The minimum principle derived herein is more general 
than both these principles, being applicable beyond the set of equilibrium positions, and also well outside the linear 
regime of non-equilibrium thermodynamics. In turn, it is based on the even broader generic formulation of the least 
action bound given herein, applicable to any system which can be analysed by Jaynes' method. 



5. CONCLUSIONS 



In this study, the manifold of stationary positions inferred by Jaynes' MaxEnt and MaxREnt principles - considered 
as a function of the moment constraints or their conjugate Lagrangian multipliers - is endowed with a Riemannian 
geometric description, based on the second differential tensor of the entropy or its Legendre transform (negative 
Massieu function) obtained from Jaynes' method. The analysis provides a generalised least action bound applicable 
to all Jaynesian systems, which provides a lower bound to the cost (in generic entropy units) of a transition between 
inferred positions along a specified path, at specified rates of change of the control parameters. The analysis therefore 
extends the concepts of "finite time thermodynamics" , developed over the past three decades [lOj [TTJ [12j [13j HH [15] 

[Miiaiisiisiiaiin 

generic Jaynes domain, providing a link between purely static (stationary) inferred positions of a system, and dynamic 
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transitions between these positions (as a function of time or some other coordinate). If the path is unspecified, the 
analysis gives an absolute lower bound for the cost of the transition, corresponding to the geodesic of the Riemannian 
hypersurface. 

The analysis is then applied to (i) an equilibrium thermodynamic system subject to mean internal energy and volume 
constraints, and (ii) a flow system at steady state, subject to constraints on the mean heat, mass and momentum fluxes 
and chemical reaction rates. The first example recovers the minimum entropy cost of a transition between equilibrium 
positions, a widely used result of finite-time thermodynamics. The second example leads to a new minimum entropy 
production principle, for the cost of a transition between steady state positions of a flow system. The analyses reveal 
the tremendous utility of Jaynes' MaxEnt and MinXEnt methods augmented by the generalised least action bound, 
for the analysis of probabilistic systems of all kinds. 
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APPENDIX A: RIEMANNIAN GEOMETRIC CONSIDERATIONS 



It is necessary to examine several salient features of the Riemannian geometric interpretation adopted herein 601 162] . 
Consider a hypersurface represented by the position vector x = [#i, x n ] T , embedded within the n-dimensional 
space defined by the coordinates (xi, ...,x n ). For analysis, this hypersurface can be converted to the parametric 
representation x(u) = [xi(u), x n (u)] T , where u = [u±, ...,^ n _i] T is the (n — l)-dimensional vector of parameters 
Uj, consisting of coordinates on the hypersurface. The first fundamental form of this geometry is defined by the metric 
EQJ|62]: 

n—l n—l 

ds 2 = dx • dx = ^ ^ aijduiduj = du T a du (Al) 

i=l j=l 

in which, by elementary calculus, the components of the tensor a can be shown to be: 

- ® X ® X ( A?} 

lJ dui duj 



Accordingly, a is symmetric. By Euclidean geometry, (Al) can be used to calculate distances between two points a 
and b on the hypersurface cc, on the path defined by u: 



L x = / ds = / V du T a du — / V du T a du d£ 

J a J a J £ a 



(A3) 



where the overdot indicates the derivative with respect to the path parameter £. The second fundamental form of the 
hypersurface is then defined by [60j [62] : 

n— 1 n—l 

— dx ■ dn = ^ ^ bijduiduj — du T b du (A4) 
i=i j=i 

where n is the unit normal vector to the hypersurface. By differential calculus, it can be shown that: 

dx 

duidu 



K = ■ 71 ( A5 ) 



The second fundamental form is not considered as a metric with which to calculate distances, but is used to examine 
the tangency and curvature properties of the manifold x [6Q| [62] . 

In the present study, we wish to adopt the Jaynesian matrix g or 7 as a Riemannian metric tensor for the calcula- 
tion of arc lengths on the .R-dimensional stationary state hypersurface, embedded in the (R + 1) -dimensional space 
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defined by (-£)*,{ (A-)}) or (0>{<V})- We therefore adopt the (somewhat peculiar) approach in which the coordinates 
[#2, #r+i] t are selected as the surface parameters [ui, , u j r] t ; i.e. with the hypersurface x%* = [$)*, (/i), (/i?)] T 
parameterised by u^* = f and with = [</>, Ai, A#] T parameterised by = A. Two necessary conditions for the 
use of g or 7 as metric tensors is that they be symmetric and positive definite (or semi-definite); since they constitute 
Hessian matrices of the concave generic entropy $)* or convex potential function 0, these conditions are satisfied, not 
only in thermodynamic applications but within the generic Jaynes formulation (with semi-definite behaviour only at 
singularities) [6j [15] . However, g and 7 are related to a second, rather than a first, fundamental form p~5j|3T]. For 
g or 7 to be considered as metric tensors, they must be able to generat e th e first fundamental form of some position 
vector which describes the hypersurface. In mathematical terms, from (Al): 



dsL = df T gdf = dufy* T a^* du%* , 



dsj 



dA T 7 dA = du<^ 3^ duA 



(A6) 
(A7) 



From ([8]), (10), (13) and (A2), taking advantage of tensor symmetries, the metric components must therefore satisfy: 

duj duj d 2 $)* 



In 



dX r 

d(fm) ' W) = d(fm)d(f r ) = d(f m ) 

dn art d 2 <\) _ d(f r ) 

dX r 



dX 



dX m dX r dX Ti 



where u;(/) and f2(A) are new .R-dimensional position vectors, which from (A10), are related by: 

a^* a^ 



I 



(A8) 
(A9) 

(A10) 



In consequence, the metrics (14)-(15) and (20)-(|2l)) and arc lengths (Jl6|-(19) used herein are not measures of distance 
on the stationary state hypersurface defined by $)*({(f r }}) or 4>({KT), but rather, on the transformed hypersurface 
given by uj or Q. In addition to the symmetry and positive definiteness conditions, it is therefore also necessary and 
sufficient that the hypersurface defined by u> or Q exists within R^, is continuous and continuously differentiate - 
at least up to first order - except in the neighbourhood of singularities. 
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